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Statement  of  the  Problem  Studied 


Advanced  structural  materials  are  generally  developed  using  various  empirical 
methods.  Such  methods  are  typically  time-consuming  and  costly  and  do  not  enable  the 
advances  in  development  of  one  type  of  material  to  be  easily  utilized  to  make  similar 
advances  in  development  of  materials  of  other  types.  The  main  objective  of  the  present 
work  was  to  develop  a  multi-length  scale  methodology  for  analyzing  the  relationships 
between:  (a)  process  parameters  and  the  resulting  materials  microstructure  and  (b) 
microstructure  properties  and  performance  in  advanced  materials  (coatings,  functionally- 
graded  materials,  etc.)  processed  by  CVD  and  LENS™  techniques.  The  work  was  also 
aimed  at  combining  various  scientific  and  engineering  disciplines  such  as:  chemical 
thermodynamics  and  kinetics,  fluid  dynamics  and  heat  transfer,  stochastic  computational 
methods  for  prediction  of  microstructural  evolution,  atomistic  simulation  techniques, 
micromechanics  approaches  and  nonlinear  fracture  mechanics. 
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Summary  of  Most  Important  Results 

1.  Atomic  simulations  combined  with  the  first-principle  electronic-structure  calculation- 
based  interatomic  potentials,  have  been  utilized  to  determine  the  grain  boundary/interface 
decohesion  potentials  and  relate  them  to  the  interface  orientation  and  grains 
misorientation.  This  approach  revealed  a  critical  role  of  grain  boundary  sliding  in  the 
process  of  decohesion.  This  work  was  done  in  collaboration  with  Dr.  Genrich  L.  Krasko  of 
the  U.S.  Army  Research  Laboratory  in  Aberdeen  Proving  Ground,  and  has  resulted  in 
eight  journal  publications  [1-8]  and  one  technical  report  [9].  Selected  results  pertaining  to 
this  portion  of  the  project  are  given  in  Appendix  A. 

2.  The  Voronoi  Cell  Finite  Element  formulation  has  been  utilized  to  determine  the  effective 
microstructure-sensitive  elastic  and  plastic  properties  of  functionally  graded  materials. 
Results  show  that  the  effective  properties  of  FGMs  are  more  reliably  determined  using  this 
method  than  the  more  commonly  used  self-consistent  approach.  This  portion  of  the  work 
has  resulted  in  two  journal  publications  [10,11].  Selected  results  pertaining  to  this  portion 
of  the  project  are  given  in  Appendix  B. 

3.  A  combined  finite  element,  singular  value  decomposition  method  has  been  established  to 
study  the  asymptotic  singular  stress  fields,  including  higher  order  terms,  for  cracks, 
interface  cracks  and  wedges  in  power  law  hardening  materials.  This  method  is  versatile 
with  respect  to  geometry  and  loading.  These  capabilities  have  led  to  the  following  two  key 
findings  for  fracture  in  functionally  graded  materials:  1)  For  plane  strain  fracture  near 
conditions  of  pure  mode  I,  the  asymptotic  stress  field  is  not  mixed  mode.  Rather  it  is  in  a 
higher  order  form.  In  FGMs,  due  to  the  combination  of  materials,  it  is  likely  that  a  crack 
will  grow  in  a  manner  that  fluctuates  around  pure  symmetric  mode  I  conditions,  where  this 
new  asymptotic  solution  is  required.  2)  For  interface  fracture  it  has  been  discovered  that 
two  asymptotic  solutions  are  necessary  to  describe  the  full  range  of  loading  conditions.  The 
first  is  a  higher  order  solution,  again  near  mode  I,  and  the  second  is  believed  to  be  non- 
separable.  Efforts  to  determine  this  non-separable  solution  are  underway.  This  portion  of 
the  work  resulted  in  five  journal  publications  [12-16]  with  one  in  review  [17],  another 
about  to  be  submitted  [18]  and  two  presentations  [19,20].  Selected  results  pertaining  to  this 
portion  of  the  project  are  given  in  Appendix  C. 

4.  Multi-length  scale  modeling  of  the  chemical  vapor  deposition  process  is  carried  out  to 
combine  the  reactor  scale,  the  grain-size  scale  and  the  atomistic  scale.  This  modeling  effort 
combines  various  approaches  from  the  disciplines  of  reactive-gas  fluid  dynamics  and  heat 
transfer,  gas  and  surface  chemical  thermodynamics  and  kinetics,  van-der  Drift-type 
modeling  of  microstructure  evolution  and  kinetic  Monte  Carlo  atomistic-scale  deposition 
modeling  techniques.  The  approach  allows  determination  of  the  relationship  between  the 
process  parameter  and  the  microstructure  and  quality  of  the  CVD-grown  film/coating. 

This  portion  of  the  work  has  resulted  in  five  journal  publications  [21-25]  and  one 
conference  presentation  [26].  Selected  results  pertaining  to  this  portion  of  the  project  are 
given  in  Appendix  D. 
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5.  To  develop  experimental  capabilities  needed  to  complement  the  work  described  above,  a 
large-scale  CVD  reactor  has  been  installed  at  Clemson  University.  This  equipment  is 
valued  at  approximately  $500,000,  and  was  donated  by  Dupont  to  the  P.I.,  who  was 
instrumental  in  the  design  of  the  equipment.  Preliminary  runs  involving  CVD  of  TiN 
coatings  have  been  carried  out  and  the  results  obtained  indicate  that  major  redesign  of  the 
CVD  reactor  is  needed  in  order  to  enable  the  level  of  process  control  suggested  by  the 
multi-length  scale  simulation  studies  after  CVD  process.  Selected  pictures  showing  the 
CVD  reactor  are  given  in  Appendix  E. 

6.  A  large-scale  parallel-computing  Monte  Carlo  Model  is  developed  to  analyze 
microstructure  evolution  during  LENS™  rapid  manufacturing  process.  The  model  enables 
establishment  of  the  relationship  between  the  process  parameters  and  the  resulting 
materials  microstructure.  This  portion  of  the  work  has  resulted  in  one  manuscript,  which 
will  be  submitted  for  publication  in  January  2001  [27]. 
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Figure  7  Equilibrium  configuration  of  the  (110) p/( 111) Y  inter¬ 
face  (bicrystal  A):  (a)  Atomic  positions  projected  onto  the 
(110)g/(211)v  plane;  (b)  Atomic  positions  projected  onto  the 
(110)p/(lll)y  interface  plane.  Interfacial  dislocations  are  marked 
as  Di  and  D2. 
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Figure  8  Equilibrium  configuration  of  the  (001)p/(011)y  inter¬ 
face  (bicrystal  B):  (a)  Atomic  positions  projected  onto  the 
(110)p/(lll)y  plane;  (b)  Atomic  positions  projected  onto  the 
(001)p/(011)y  interface  plane.  Interfacial  dislocations  are  marked 
as  Dx  and  D2. 
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Figure  9  Equilibrium  configuration  of  the  (110)p/(211)v  interface 
(bicrystal  C):  (a)  Atomic  positions  projected  onto  the  (110)p/(lll)Y 
plane;  (b)  Atomic  positions  projected  onto  the  (001)p/(011)Y  plane; 
(c)  Atomic  positions  projected  onto  the  (110)p/(211)Y  interface 
plane.  Interfacial  dislocations  are  marked  as  Di  and  D2. 
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Figure  10  Equilibrium  configuration  of  the  (lll)yi/(lll)y2  inter¬ 
face  (bicrystal  D):  (a)  Atomic  positions  projected  onto  the 
(011)vl/(011)v2  plane;  (b)  Atomic  positions  projected  onto  the 
(lll)yi/(lll)Y2  interface  plane.  Interfacial  screw  dislocations  are 
marked  as  Di. 
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Figure  11  Variation  of  the  decohesion  potential  for  the  (H0)r/(111)y 
interface  with  the  normal  displacement  Un  in  the  [H0]p/[lTI]Y  direc¬ 
tion  and  tangential  displacements  Ut  and  Ub  in  the  [001]p/[0ll]Y  and 
[110]p/[211]Y  directions,  respectively. 
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Figure  12  Variation  of  the  decohesion  potential  for  the  interface 
(001)p/(011)y  with  the  normal  displacement  Uu  in  the  [001]e/[011]v 
direction  and_tangential  displacements  Ut  and  Ub  in  the  [110]p/[lll]Y 
and  [110]p/[211]y  directions,  respectively. 


19 


Figure  13  Variation  of  the  decohesion  potential  for  the  interface 
(110)p/(211)Y  with  the  normal  displacement  Un  in  the  [ll0]p/[2lI]Y 
direction  and  tangential  displacements  Ut  and  Ub  in  the 
[110]p/[lll]Y  and  [001]p/[011]Y  directions,  respectively. 
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Figure  14  Variation  of  the  decohesion  potential  for  the  interface 
(lll)vl/(lll)v2  with  the  normal  displacement  Un  in  the 
[111]yi/[111]v2  direction  ajod  tangential  displacements  Ut  and  Ub  in 
the  [011]yi/[011]72  and  [211]Yi/[211]Y2  directions,  respectively. 
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Table  I:  Parameters  Characterizing  Interfacial  Dislocation  for  the  three  Beta/ 
Gamma  and  one  Gamma/Gamma  Interfaces  Analyzed  in  the  Present  Work 


(110)p/(lll)y  Interface:  Bicrystal  A 

Parameter 

Dislocation  Type 

Di 

Line  Direction, 

1 

[l  T  o]p 

[2  ll]y 

[0  o  i]p 

[0  1  l]y 

Burger’s  Vector, 

b 

1/2  aB  [0  0  l]p 

1/2  ayCyA/  ay2+Cy2  [0  1  l]y 

4  aBV2  [T 1  0]B 

4  ayCy/v/ayZ+5Cy2  [2  1  l]y 

Dislocation  Spacing, 

X 

7/2  ag 

4  fly  Os// \/ 

lOag/y/2 

14&yCv/ v/  S^/^+SCv" 

i  1  v  i  J 

(001)  p/(01I)y  Interface:  Bicrystal  B 

Parameter 

Dislocation  Type 

Di 

D2 

Line  Direction, 

1 

[1 1  0]p 
[1  1  1]Y 

[3  4  0]p 
[3  2  2]y 

Burger’s  Vector, 

b 

2a6/v/2[TlO]B 

2  ayCy/^/a^^+SCy 2  [2  1  l]y 

aB/5[4  3  0]B 

ayC,,/v/9ayZ+25cyZ  [4  3  3],, 

Dislocation  Spacing, 

X 

7aBV2 

10  SkyCyfy/ fly'^+SCy'^ 

21  aB/5 

20  ayCy/v7  9a,/+25cyZ 

(110)p/(211)y  Interface:  Bicrystal  C 

Parameter 

Dislocation  Type 

Di 

I>2 

Line  Direction, 

1 

[1  1  0]p 
[1  1  1]Y 

[3  3  7]  p 
[1  1  3]y 

Burger’s  Vector, 
b 

1/2  ap  [0  0  l]p 

1/2  ayCy/-v/ayZ+CyZ  [0  1  l]y 

aBVT34[7  7  6]B 

3yC  yj\/  3y^H"65Cy^  [471]y 

Dislocation  Spacing, 

X 

7/2  aB 

4  ftyCy f\/  3y^+Cy^ 

21  aB/v/134 

22  3yCy / y/  fly^+^SCy^ 

(110)p/(211)y  Interface:  Bicrystal  D 

Parameter 

Dislocation  Type 

Di 

Line  Direction, 

1 

[0Tl]y 

Burger’s  Vector, 
b 

1/2  ayCy/v/ayZ+Cyi  [0  1  l]y 

Dislocation  Spacing, 

X 

25/2  ay  v"  (ayz+2cyZ)/(aYz+CyZ) 
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Table  II.  Decohesion  Potential  Parameters  for  the  Three  Beta/Gamma  and 
One  Gamma/Gamma  Interface  Analyzed  in  the  Present  Work 


Parameter 

Interface  Type 

(U0)p/(111)Y 

(001)  p/ (Oil)  y 

(110)p/(211)y 

(Hl)Yl/(lll)Y2 

Designation 

A 

B 

C 

D 

n-direction 

[Ho]p/[m]Y 

■eiatsaaaf 

IBBiBIB 

t-direction 

[HO]  p/[lll]  y 

UnWlB 

b-direction 

6n>  A 

0.5 

0.5 

0.5 

0.5 

3.21 

4.54 

6.81 

123.5 

A 

4.54 

4.54 

3.21 

5.78 

a0 

-0.013 

-0.21 

-0.148 

-0.039 

ai 

-1.5 

-1 

0.42 

-1 

o-i 

0 

0 

-0.2 

0 

Ct3 

0 

0 

-3.42 

0 

a4 

0.5 

0 

1.17 

0 

05 

0 

0 

-5 

0 

^max?  GPa 

7.37 

4.57 

7.29 

2.01 

^max,t?  GPa 

1.09 

11.77 

29.2 

0.160 

1.11 

5.60 

1.20 

1.65 

<t>(Un-+oo)5  J/m2 

2.02 

1.25 

2.03 

0.561 
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Figure  9  Contour  plots  at  the  overall  normal  strain  of  2%  in  the  vertical  direction  of:  (a)  the  normal 
equivalent  plastic  strain  (in  percents),  (b)  the  changes  in  Euler  angle  (in  degrees)  and  (c)  the  hydro¬ 
static  stress  (in  lOOMPa)  for  the  single-phase  (matrix)  polycrystalline  material  in  which  the  grain 
boundaries  are  represented  using  the  cohesive  zone  model,  and  (d),  (e)  and  (f)  are  the  corresponding 
contour  plots  for  the  two-phase  (matrix  +  stable  beta)  material. 
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Figure  10  Contour  plots  for  the  normal  equivalent  plastic  strain  (in  percents)  for  the  single-phase 
(matrix)  material  at  the  overall  normal  strains  of  (a)  1  %,  (b)  1.7%  and  (c)  2.1%.  The  corresponding 
contour  plots  for  the  change  in  Euler  angle  (in  degrees)  (d),  (e)  and  (f).  The  corresponding  contour 
plots  for  the  hydrostatic  stress  (in  lOOMPa),  (g),  (h)  and  (i). 
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Appendix  B 


Voronoi  Cell  Finite  Element  Analysis 
Of 

Effective  Properties  of  Functionally  Graded  Materials 
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(b)  (c)  (d) 


Figure  5.  (a)  Schematic  representation  of  the  microstructure  of  a  graded  ceramic/ 
metal  region.  The  microstructure  consisting  of  dispersed  particles  embedded  in 
a  continuous  matrix  is  observed  for  large  and  small  volume  fractions  of  the  metal 
and  ceramic,  as  in  (b)  and  (d).  In  the  intermediate  region  where  the  volume  frac¬ 
tions  of  the  two  materials  are  comparable  to  each  other,  the  microstructure  con¬ 
sists  of  intertwined  clusters  of  the  two  phases. 
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Figure  8.  Effect  of  the  number  of  Voronoi  cell  elements  on  the  range  of  in-plane  prop¬ 
erties  for  the  microstructure  consisting  of  0.2  volume  fraction  of  elliptical  inclusions 
(aspect  ratio  in  the  range  1-2).  One  of  the  100  element  square  RMEs  used  is  shown 
in  (a) ,  while  one  of  the  25  element  square  RMEs  used  is  shown  in  Figures  5  (a)  and  5  (b) . 


100 

Number  of  Elements 


Figure  9.  Effect  of  the  number  of  Voronoi  cell  elements  on  the  range  of  in-plane  prop¬ 
erties  for  the  microstructure  consisting  of  intertwined  clusters  of  the  two  materials 
(volume  fraction  of  the  two  materials  —0.5).  One  of  the  600  element  RMEs  used  is 
shown  in  (a). 
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MgO  Volume  Fraction  MgO  Volume  Fraction 


Figure  10.  Comparison  of  the  experimentally  measured  Young’s  modulus  (a)  and  Pois¬ 
son’s  ratio  (b)  inNi/MgO  FGM  with  the  ones  predicted  by  the  Self  Consistent  Method 
(SCM),  the  Equivalent  Inclusion  Method  (EIM)  and  the  Yoronoi  Cell  Finite  Element 
Method  (VCFEM).  The  experimentally  measured  porosities  are  shown  in  the  insert 
of  figure  (b). 


TIC  Volume  Fraction 


TiC  Volume  Fraction 


Figure  1 1 .  Comparison  of  the  experimentally  measured  Young’s  modulus  (a)  and  Pois¬ 
son’s  ratio  (b)  in  N^Al/TiC  FGM  with  die  ones  predicted  by  the  Self  Consistent  Meth¬ 
od  (SCM),  the  Equivalent  Inclusion  Method  (EIM)  and  the  Voronoi  Cell  Finite  Ele¬ 
ment  Method  (YCFEM).  The  experimentally  measured  porosities  are  shown  in  the 
insert  of  figure  (b). 


Figure  10.  The  finite  element  simulation  results  of  thermally  induced:  (a)  interface  decohesion,  (b) 
cracks  in  AI2O3  (layer  1),  90%Al203/10%  316  stainless  steel  (layer  2)  and  70%Al2C>3/30%  316 

I 

stainless  steel  (layer  3),  (c)  plastic  strain  and  (d)  voids  in  316  stainless  steel,  (e)  radial  residual  stress, 
Oj,  and  (f)  axial  residual  stress,  ct2,  in  the  case  of  the  graded  316  stainless  steel/Al203  interface  char¬ 
acterized  by  the  exponent  p=2.0. 1/2, 2/3,  etc.  in  (a)  refer  to  the  interfaces  separating  layers  1  and 
2, 2  and  3,  etc. 
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Appendix  C 


Multiple  Asymptotic  Solutions 
For  Mixed  Mode  Fracture 
In  Power  Law  Hardening  Materials 
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This  figure  shows  how  the  linear  elastic  double  root  for  plane  strain 
fracture  in  the  upper  figure  (a  =  360° ,  Sj  =  s2  =  -0.5 )  splits  into  two 
distinct  asymptotic  solutions  when  material  hardening  is  taken  into 
account  in  the  lower  figure.  The  values  of  si  and  S2  are  defined  by, 

CTjj  (r,  0)  =  Ajrs‘  si  (0)  +  A2rSj  a?  (0)  .  The  mode  I  dominant  solution  is 

higher  order  (s2  >  s, ),  while  the  mode  II  dominant  solution  is  mixed 
mode(s2=s,). 


36 


-O.Oi 


Figure  7a. 

n=l,  linear  elastic, 
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This  figure  shows  how  the  linear  elastic  double  root  for  plane  strain  interface  fracture 
in  the  upper  figure  (a  =  l  80°  ,Sl=S2=  -0.5 )  splits  into  two  distinct  asymptotic 
solutions  when  material  hardening  is  taken  into  account  in  the  lower  figure.  The 
values  of  S\  and  s2  are  defined  by,  c;j  (r,  0)  =  1 s“  a?.  (0)  +  A2r 52  a?  (0) .  The  mode  I 

dominant  solution  is  higher  order  (  s2  >  sj,  while  the  mode  II  dominant  solution  is 
believed  to  be  non-separable. 
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Multi-length  Scale  Modeling 
0f 

Chemical  Vapor  Deposition  Process 
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Figure  2.  Sequence  of  steps  involved  in  the  nucleation  of  a  new  layer  on  an 
atomically  flat  (111)  surface.  Open  circles  represent  diamond 
carbon  atoms,  shaded  circles  designate  hydrocarbon  carbon  atoms, 
and  small  black  circles  stand  for  hydrogen  atoms.  Surface 
reactions  (listed  in  Table  I)  involved  in  various  steps  of  the  process 
are  given  in  parenthesis. 
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Reactions  Involved 
2x[(3)  or  (-7)] 


Growth  of  an  (lOO)-oriented  film  by  trough  insertion  mechanism. 
Open  circles  represent  diamond  carbon  atoms,  shaded  circles 
designate  hydrocarbon  carbon  atoms,  and  small  black  circles  stand 
for  hydrogen  atoms.  Surface  reactions  (listed  in  Table  I)  involved 
in  various  steps  of  the  process  are  given  in  parenthesis. 


Deposition  Rate,  (um/hr) 


Figure  5 .  Reactor-  and  atomic  scale  analyses  predicted  (111)-  and  ( 1 00)- 
oriented  film  deposition  rates  as  a  function  of  (a)  the  substrate 
temperature  and  (b)  concentration  of  CH4  in  feed  gas.  The 
remaining  processing  condition  are  as  indicated  in  Figure  4. 
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Figure  5. 


Top  view  of  an  (lll)-onented  diamond  film  obtained  under 
the  following  CVD  conditions  in  the  reactor:  Reactive  gas  at 
the  reactor  inlet  (0.4%  CH4,  92.5%  H2),  Theator=  2000  K, 
Tsubstrate=  1000  K,  p=  20.25  Ton;  Heater-to-Substrate 
Distance  =  1.3cm.  Deposition  times:  (a)  0.87s;  (b)  1.81s;  (c) 
2.07s  and  (d)  2.85s.  Nomenclature:  B  -  3-carbon  bridge, 

C  -  Twin  covered  by  regular  crystal,  D  -  Dislocation  loop, 
E  -  Edge,  G  -  Gap,  I  -  Island,  K  -  Kink,  N  -  Nucleus, 

T  -  Twin,  V- Void 


1 


Figure  6.  The  side  view  of  the  four  (11  l)-oriented  diamond  films 
shown  in  Figures  5(a)-(d). 


43 


Figure  9.  Top  view  of  four  (100)  surface  configurations  obtained  under 
the  CVD  conditions  identical  to  the  ones  listed  in  Figure  5. 
Deposition  times:  (a)  0.01s;  (b)  0.018s;  (c)  0.032s  and  (d) 
0.208s.  Nomenclature:  B — BCN  mechanism,  D — Dimer 
insertion  mechanism,  P — Pit,  I — Island,  T — Trough 
insertion  mechanism. 
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Figure  10.  Side  view  of  the  four  (lOO)-oriented  diamond  films  shown  in 
Figures  9(a)-(d). 


45 


Concentration  of  Point  Defects,  (at%) 


Figure  8.  Effect  of  the  substrate  temperature  on  the  concentration  of 
vacancies  and  entrapped  hydrogen  atoms  in  (111)-  and 
(lOO)-oriented  diamond  films,  (a),  and  twins  in  (lll)-oriented 
diamond  films,  (b),  under  the  CVD  conditions  specified  in 
Figure  5.  Error  bars  represent  on  standard  deviation  over  five 
simulation  runs. 
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Normalized  Film  Thickness  Normalized  Film  Thickness 
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The  microstructure  of  polycrystalline  diamond  films  deposited  under  the 
growth-rate  parameter,  a=:  (a)  1.0;  (b)  1.05;  (c)  1.5;  (d)  2.95  and  (e)  3.0. 
Both  the  x-  and  y-axis  are  normalized  with  respect  to  the  average 
nucleus  spacing,  d©.  The  following  line  type  nomenclature  is  used:  dash- 
dot=grain  boundaries,  solid={100}  facets,  dotted  ={111}  facets,  dashed  = 
{100}/{100}  facets,  and  three-dot  space={lll}/{lll}  facets. 


Figure  4.  (a)  The  distribution  function  for  the  grain  size  normalized 

with  respect  to  the  average  nucleus  spacing,  do,  at  five  film 
thicknesses  for  the  growth-rate  parameter  a=1.5.  (b)  The 
distribution  function  for  the  grain  size  normalized  with 
respect  to  the  corresponding  average  grain  size,  d,  for  five  a 
values  and  three  film  thicknesses. 
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Appendix  E 


Bernes  2000 
Prototype  Scale 

Chemical  Vapor  Deposition  Reactor 
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